home *** CD-ROM | disk | FTP | other *** search
/ FishMarket 1.0 / FishMarket v1.0.iso / fishies / 051-075 / disk_065 / prep / sieve.f < prev    next >
Text File  |  1992-05-06  |  558b  |  25 lines

  1.       integer f(0:8190 )
  2.       integer i, p, k, c, n
  3.       do 12500 n = 1, 10
  4.          c = 0
  5.       do 10000 i000 = (  0), ( 8190 ), 1
  6.          f(i000) = 1
  7. 10000 continue
  8.       do 12501 i = 0, 8190
  9.             if ( f(i) .ne.  0 ) then
  10.                p = 2*i + 3
  11.                k = i + p
  12. 17500 continue
  13.       if (.not.( k .le.  8190  )) goto 15000
  14.                   f(k) = 0
  15.                   k = k + p
  16.       goto 17500
  17. 15000 continue
  18.                c = c + 1
  19.             end if
  20. 12501 continue
  21. 12500 continue
  22.       write(*,*) c, ' primes'
  23.       stop
  24.       end
  25.